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Recent developments in string theory have led to 5-dimensional warped spacetime models in 
which standard- model fields are confined to a 3-brane (the observed universe), while gravity can 
propagate in the fifth dimension. Gravity is localized near the brane at low energies, even if the extra 
dimension is noncompact. A review is given of the classical geometry and dynamics of these brane- 
world models. The field equations on the brane modify the general relativity equations in two ways: 
local 5-D effects are imprinted on the brane as a result of its embedding, and are significant at high 
energies; nonlocal effects arise from the 5-D Weyl tensor. The Weyl tensor transmits tidal (Coulomb) , 
gravitomagnetic and gravitational wave effects to the brane from the 5-D nonlocal gravitational field. 
Local high-energy effects modify the dynamics of inflation, and increase the amplitude of scalar and 
tensor perturbations generated by inflation. Nonlocal effects introduce new features in cosmological 
perturbations. They induce a non-adiabatic mode in scalar perturbations and massive modes in 
' vector and tensor perturbations, and they can support vector perturbations even in the absence 

of matter vorticity. In astrophysics, local and nonlocal effects introduce fundamental changes to 
gravitational collapse and black hole solutions. 

t"*- 1 I. INTRODUCTION 

7— i ; 

\ At high enough energies, Einstein's theory of general relativity breaks down and is likely to be a limit of a more 
■ general theory. In string theory/ M theory, gravity is a truly higher-dimensional theory, becoming effectively 4- 
0^ ' dimensional at lower energies. Recent developments may offer a promising road towards a quantum gravity theory [[lj . 
In brane-world models inspired by string/M theory, the standard-model fields are confined to a 3-branc, while the 
gravitational field can propagate in 3 + d dimensions (the 'bulk'). The d extra dimensions need not all be small, or 
even compact: recently Randall and Sundrum Q have shown that for d = 1, gravity can be localized on a single 
3-brane even when the fifth dimension is infinite. This noncompact localization arises via the exponential 'warp' 
factor in the non-factorizable metric: 



o 
o 



o 
o 

^ dS 2 = exp(-2|y|/^) [-dt 2 + dx 2 ] + dy 2 . (1) 

For y 7^ 0, this metric satisfies the 5-dimensional Einstein equations with negative 5-dimensional cosmological constant, 
J^P| A cx — 1~ 2 . The brane is located at y = 0, and the induced metric on the brane is a Minkowski metric. The bulk is a 
^ i 5-dimensional anti-de Sitter metric, with y = as boundary, so that y < is identified with y > 0, reflecting the 
I^j ' symmetry, with the brane as fixed point, that arises in string theory. 

Perturbation of the metric ([!]) shows that the Newtonian gravitational potential on the brane is recovered at lowest 
' order: 

„, , GM ( 2£ 2 \ 

Thus 4-dimensional gravity is recovered at low energies, with a first-order correction that is constrained by current 
sub- millimetre experiments Q . The lowest order term corresponds to the massless graviton mode, bound to the brane, 
while the corrections arise from massive Kaluza-Klein modes in the bulk. Generalizing the Randall-Sundrum model 
to allow for matter on the brane leads to a generalization of the metric ([!]), and to a breaking of conformal flatness, 
since matter on the brane in general induces Weyl curvature in the bulk. Indeed, the massive Kaluza-Klein modes 
that produce the corrective terms in Eq. (|J) reflect the bulk Weyl curvature that arises from a matter source on the 
brane. 

At a classical level, the brane-world models have a rich geometrical structure, in which the bulk Weyl curvature 
tensor and its interaction with the curvature of the brane play an important role. The classical geometry and dynamics 
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of brane-world models are important because they provide a foundation for testing the astrophysical and cosmological 
predictions of these models against observations. Any modification of Einstein's theory that arises in the quest for 
quantum gravity needs to pass the tests provided by increasingly accurate observations. Observers are necessarily 
bound to the brane, unable to access the bulk. Bulk effects are felt indirectly via their imprint on the brane geometry 
and dynamics. A covariant Lagrangian approach to the brane-world M brings out clearly the role of the Weyl tensors, 
the kinematical and dynamical quantities, and their relation to observations. 



II. FIELD EQUATIONS ON THE BRANE 



Instead of starting from a metric ansatz in special coordinates, one can generalize the Randall-Sundrum model by a 
covariant geometric approach, given by Shiromizu, Maeda and Sasaki |J. (See also || for related geometrical results). 
The unit normal to the brane n A defines the induced metric on the brane (and on all hypersurfaces orthogonal to 

n A ), 

9ab = 9ab - n A n B , (3) 

where we use tildes to denote the 5-dimensional generalization of standard general relativity quantities. Without loss 
of generality, we can take n A to be geodesic. A natural choice of coordinates is x A — (2^,2/), where x 11 = {t,x z ) are 
spacetime coordinates on the brane and ua — Sa v ■ The extrinsic curvature orthogonal to n A is 

K A b = \Cn9AB = gA C ^criB , (4) 

so that K[ab] — = Kab^ b , where square brackets denote anti-symmetrization and C is the Lie derivative. 
The Gauss equation gives the 4-dimensional curvature tensor as 

Rabcd — Refgh 9a 9b 9C9D H + 2K A[C K D]B , (5) 
and the Codazzi (or Mainardi-Codazzi) equation determines the change of Kab '■ 

V B K B A - V A K B B = Rbc 9A B n c . (6) 
The 5-dimensional Einstein equations are 



Gab = k 



^gAB + <%) {-XgAB + Tab} , (J) 



where k 2 = 87r/Mp, with M p the fundamental 5-dimensional Planck mass, which is typically much less than the 
effective Planck mass on the brane, M p = 1.2 X 10 19 GeV. The brane tension is A, and fields confined to the brane 
make up the brane energy-momentum tensor Tab, with Tab^ 3 — 0. The Randall-Sundrum model may be further 
generalized by allowing for a scalar field in the bulk |jj . 
Using Eqs. (|) and (@), it follows that §] 

Gab = ~\T?lg A B + K C C K AB ~ K A C K CB 

+ \\k cd K cd -(K c c ) 2 ]9ab~£ab, (8) 



where 



Sab = Cacbd n c n D , (9) 



is the projection of the bulk Weyl tensor orthogonal to n A , with £[ab] = = Ea A ■ Evaluating Eq. (|8|) on the brane 
(strictly, as y — > ±0) will give the field equations. First, we need to determine Kab at the brane. The junction 
conditions across the brane imply that gAB is continuous, while Kab undergoes a jump due to the energy-momentum 
on the brane: 

K\ B - K AB = -k 2 [Tab + |( A - T c°) 9ab] ■ (10) 

The Z2 symmetry implies that 
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K AB = -K+ B , (11) 

and then 

K AB = [T AB + M A ~ T c C ) 9ab] , (12) 

where we have dropped the (+) and we evaluate quantities on the brane by taking the limit y — > +0. 
Finally we arrive at the induced field equations on the brane Q : 

Gfj, u — —Ag^v + K^T^y + k S^ v — E^ v , (13) 

where k 2 = 8ir/M 2 . The energy scales are related to each other via 

A = 6^ , A = \H 2 (X + i« 2 A 2 ) . (14) 

The higher-dimensional modifications of the standard Einstein equations on the brane are of two forms: first, the 
matter fields contribute local quadratic energy- momentum corrections via the tensor S^ v , which arise from the extrinsic 
curvature, and second, there are nonlocal effects from the free gravitational field in the bulk, transmitted via the 
projection £ „„ of the bulk Weyl tensor. The local corrections are given by 

S [iu 12^°^ ^fJ-^ 4-^/iaT v 24^M^ ^TocpT ^ (T a ) . (1^) 

The Weyl tensor Cabcd represents the free, nonlocal gravitational field in the bulk. The local part of the bulk 
gravitational field is the Einstein tensor Gab, which is determined locally via the bulk field equations (Q). Thus 
£ M „ transmits nonlocal gravitational degrees of freedom from the bulk to the brane, including tidal ( or Coulomb ), 
gravito-magnetic and transverse traceless (gravitational wave) effects. 

There may be other branes in the bulk. Branes interact gravitationally via any Weyl curvature that they generate. 
On the observer's brane at y = 0, the presence of other branes is felt indirectly through their contribution to £^ v . 

As a consequence of the Codazzi equation ([}]), the form of the bulk energy-momentum tensor (which means that 
Rbc gA B n = 0) and Z2 symmetry, it follows that the brane energy-momentum tensor is conserved: 

V"T M!/ = 0. (16) 

When there are scalar and other fields in the bulk, this is no longer in general true @,||, and non-conservation of 
T^ v reflects an exchange of energy-momentum between the brane and the bulk. In the case when there is only a 
cosmological constant in the bulk, there is no such exchange. Using Eq. ( [l6| ) in Eq. (|l^), the contracted Bianchi 
identities on the brane, V*G M „ = 0, imply that the projected Weyl tensor obeys the constraint 



= ^- WS^ . (17) 



This shows that £ M „ is sourced by energy-momentum terms, which in general include spatial gradients and time 
derivatives. Thus evolution and inhomogeneity in the matter fields can generate nonlocal gravitational effects in the 
bulk, which then 'backreact' on the brane. 

The dynamical equations on the brane are equations (|l^), ([l6]) and (|l7]). It is important to note that these equations 
are not in general closed on the brane ||, since Eq. ( fl7| ) does not determine £ M „ in general, as further discussed below. 
This reflects the fact that there are bulk degrees of freedom which cannot be predicted from data available on the 
brane, for example, incoming gravitational radiation which impinges on the brane. One needs to solve the field 
equations in the bulk in order to fully determine £^ v on the brane. 



III. COVARIANT INTERPRETATION OF GRAVITY LOCALIZATION 



In the Randall-Sundrum model, with a flat brane, localization of gravity at the brane is understood perturbatively 
via Eq. (^|). For a general matter distribution on a curved brane, we can provide a qualitative non-perturbative and 
covariant interpretation of gravity localization via tidal acceleration. 

Consider a field of observers on the brane (e.g., observers comoving with matter) with 4- velocity u^, and let u A be 
an extension off the brane (the result does not depend on the extension), so that u a ua = 0, u a ua = —1. The tidal 
acceleration in the n A direction measured by the observers is —nAR A BCDU B n c u D . Now 



3 



Rabcd — Cabcd + § \ 9a\cRd\b + 9b[dRc\a\ — \R9a[c9d]b , (18) 
so that by the field equation (^) (and recalling that Tab^ 3 = 0), 

- R A BCDn A u B n c u D = -£ AB u A u B + ±k 2 A . (19) 
Taking the limit y — > +0, we get [|j 

tidal acceleration in off-brane direction = ^k 2 A — E^u^u 1 ' . (20) 

Since A < 0, it contributes to acceleration towards the brane. This reflects the confining role of the negative bulk 
cosmological constant on the gravitational field in the generalized Randall-Sundrum type models. Equation (|2(]) also 
shows that localization of the gravitational field near the brane is enhanced if E ilv u il u v > 0, which corresponds to a 
negative effective energy density on the brane from nonlocal bulk effects. 



IV. COVARIANT DECOMPOSITION OF LOCAL AND NONLOCAL BULK EFFECTS 

The general form of the brane energy-momentum tensor for any matter fields (scalar fields, perfect fluids, kinetic 
gases, dissipative fluids, etc.), including a combination of different fields, can be covariantly given as 

T^y = pu^Uy + phfiu + 7r M „ + q^Uu + q v u^ . (21) 

Here p and p are the energy density and isotropic pressure, and — + u^u u projects orthogonal to on the 
brane. The energy flux obeys = q/^\ , and the anisotropic stress obeys -k^ v — "Kt^v) , where angled brackets denote 
the projected, symmetric and tracefree part: 

V w = h/V v , W {fiu) = [h {fl a h v f - ±/i Q/ V] W af} , (22) 

with round brackets denoting symmetrization. In an inertial frame at any point on the brane, we have — S^o and 
/i M „ = diag(0, 1, 1, 1), = (0, tt^q = 0. 

The tensor S^ u , which carries local bulk effects onto the brane, may then be irreducibly decomposed as 

S„v = i [2p 2 - 3n a ^] Ull u v + i [2p 2 + App + ir af3 TT a P - Aq a q a ] 

~ 2 P) 7r ^ + Mm 7 ^}" + Utilu) + \p%^u v) - j2q a TT a ^u u) . (23) 

This simplifies for a perfect fluid or minimally-coupled scalar field: 

S^v = j2P 2 u fl u l , + ±p (p + 2p) . (24) 

The quadratic energy-momentum corrections to standard general relativity will thus be significant for ~K A p 2 > 12k 2 /5, 
i.e., in the high-energy regime 

'> A ~(§)V < 25 > 

The lower bound arising from current tests for deviations from Newton's law is JoJ 

M p > 10 5 TcV , A 1/4 > 100 GeV . (26) 

(A much weaker limit is imposed by nucleosynthesis constraints.) 

Nonlocal effects from the bulk are encoded in the brane tensor £^ u , which can be decomposed as 

£p.v = — %r [W (u^u u + \h^ v ) + + Q^u u + Q„w M ] . (27) 

The factor 6/k 2 A, which is (k/k) 4 by Eq. (|l4|), is introduced for dimensional reasons, and also since it ensures that 
in the general relativity limit, A -1 — > 0, we have £ M „ — > 0. We have written £ M „ as an effective energy-momentum 
tensor: the bulk Weyl tensor imprints on the brane an effective energy density, stresses and energy flux. 
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The effective noniocal energy density on the brane, arising from the free gravitational field in the bulk, is 

U = -{\k 2 \)£^u»u u . 

This nonlocal energy density need not be positive [p~0| p~2[ . Since 8^ v is tracefree, the effective nonlocal pressure is 
^U. There is an effective nonlocal anisotropic stress 

on the brane, arising from the free gravitational field in the bulk, and 

is an effective nonlocal energy flux on the brane, arising from the free gravitational field in the bulk. 

If the bulk is anti-de Sitter (AdSs), as in the Randall- Sundrum model, then £^ v = 0, since the bulk is conformally 
flat: 

AdS 5 bulk: 8^=0. (28) 

The Randall-Sundrum model has a Minkowski brane, but AdSs can also admit a Friedmann brane, and satisfy the 
Einstein equations (Q). However, the most general solution with a Friedmann brane is Schwarzschild-anti de Sitter 
spacetime Jl3[ . Then it follows from the Friedmann symmetries that 

SAdS 5 bulk, Friedmann brane: Q M = = 7> M „ , (29) 

where U — only if the mass of the black hole in the bulk is zero. The presence of the black hole leads to a 'dark 
radiation' term in the Friedmann equation (see below). For a static spherically symmetric brane (e.g. a static stellar 
interior or the exterior of a nonrotating black hole) |ll|] : 

static spherical brane: Q M = . (30) 

This condition also holds for a Bianchi I brane [Q. 

The local and nonlocal bulk corrections may be consolidated into an effective total energy density, pressure, 
anisotropic stress and energy flux, since the modified field equations ( |l3| ) take the standard Einstein form with a 
redefined energy-momentum tensor: 



where 



G^v = -A 5|t „ + k 2 T^ , (31) 



tot 6 1 

T]iu H r*-V i/ ' (32) 



Then it follows from Eqs. (§|) and © that 



<° t0t = P + A ( 2 P 2 - + ~T^ U ( 33 ) 



4A " ^ > ' k 4 A 

p tot =P+j^ (V + *PP + *n»*" v - V) + ^ (34) 



~tot 

7T„„ = 7T 



1 r , . . v „ . n 6 



~{P + 3p)7r MV + TT a ^TT v ) a + q{^q v )\ + -JT^V (35) 



C = 9m + ^ (4Wm - Tr^g") + ^ • (36) 

These general expressions make clear the local and nonlocal effects. They simplify in the case of a perfect fluid (or 
minimally coupled scalar field, or isotropic one-particle distribution function), i.e., for = = n^. However, we 
note that the total energy flux and anisotropic stress do not vanish in this case in general: 

Nonlocal bulk effects can contribute to effective imperfect fluid terms even when the matter on the brane has perfect 
fluid form. 
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V. LOCAL AND NONLOCAL CONSERVATION EQUATIONS 



The brane energy-momentum tensor and the consolidated effective energy-momentum tensor are both conserved 
separately, by virtue of Eqs. (|l6|) and (|l7|). Conservation of gives the standard general relativity energy and 
momentum conservation equations 

p + Q(p + p) + D% + 2A% + = , (38) 

5( M > + 5©?/, + d mP + + P) A n + D" 71 "^ + 

+ a tu ,^-[w,q] M =0. (39) 

In these equations, an overdot denotes u"V„, O = V M u M is the volume expansion rate of the u M congruence, — 
= ^4( M ) is its 4-acceleration, g^ v = D/^it„\ is its shear rate, and ui^ — — ^curlu^, = cjaa is its vorticity rate. On a 
Friedmann brane, — lu^ = = and = 3H, where H = a/a is the Hubble rate. 
The covariant spatial curl is given by | Q 

curlF M = e^D a F /3 , curl = e aP(tl D a W p v) , (40) 

where £ M „o- is the projection orthogonal to of the brane alternating tensor, and is the projected part of the 
brane covariant derivative, defined by 

DpF '" ...p = (V /1 F a "\.. /9 )_ L - h fl u h a 1 ■ ■ ■ hffV v ri"\..s . (41) 

The covariant cross-product is [V, Y]^ — s^ a pV a Y^ . In a local inertial frame at a point on the brane, with = S^o, 
we have: = A = u; = a ^ = e ^ = [V, Y] = curl V = curl W 0pi and D^F" ".. t/ g = 5 li i 5 a j ■ ■ ■ Sp k S/iF j - ... k . 

The conservation of T*°* gives, upon using Eqs. (p3[)— (39), what we can call nonlocal conservation equations [Q. 



The nonlocal energy conservation equation is a propagation equation for U : 
u + leu + D"^ + 2^Q M + cr^TV 

= ^K 4 [6^77^ + 6(p + pJtT^TT^ + 26 (2<z% + 7r"%„) + 2A W 

- 4<f D M p + ^D%„ + 7^"D^ - 2^^^° - 2^%^] . (42) 
The nonlocal momentum conservation equation is a propagation equation for Q^: 

Q(n) + f ©Q M + §D P M + lUA^ + WV^ + A V V MV + GfivQ" - [w, Q] M 
= 2^k 4 [-4(p + p)D M p + 6(p + p)D%„ + g"7r <M1/) + 7r/D„(2p + 5p) 
- f tt"' 3 (D M 7r a/3 + 3D Q ^ M ) - 37r„ a B n af} + fq u D„q„ 
+ ipA"n^ - Sn^Apir^ + %A^ a ^ a p - n^a^qp 
+ o-^a^qp + T^fiv [w,q] v - e^ap^Tr^q,, + 4(p + p)9g M 

+ 6<^% + ^A^fr + A qi/ a a ^ a p] . (43) 



All of the matter source terms on the right of these two equations, except for the first term on the right of Eq. MS 
are imperfect fluid terms, and most of these terms are quadratic in the imperfect quantities q^ and ir^. For perfect 
fluid matter, only the D M p term on the right of Eq. ( [43] ) survives, but in realistic cosmological and astrophysical 
models, further terms will survive. For example, terms linear in tt^ will carry the photon quadrupole in cosmology 
or the shear viscous stress in stellar models. 

In general, the 4 independent equations determine 4 of the 9 independent components of E^v on the brane. What is 
missing, is an evolution equation for V^v (which has up to 5 independent components). Thus in general, the projection 
of the 5-dimensional field equations onto the brane, together with matching, does not lead to a closed system. Nor 
could we expect this to be the case, since there are bulk degrees of freedom whose impact on the brane cannot be 
predicted by brane observers. Our decomposition of £^ t „ has shown that the evolution of the nonlocal energy density 
and flux ( carrying scalar and vector modes from bulk gravitons ) is determined on the brane, while the evolution of the 
nonlocal anisotropic stress (carrying tensor, as well as scalar and vector, modes) is not. 

In special cases the missing equation does not matter. For example, if V^v = 0, as in the case of a Friedmann 
brane, then the evolution of is determined by Eqs. d42] ) and (f43|). If the brane is stationary (with Killing vector 
parallel to u M ), then evolution equations are not needed for £y, v . However, small perturbations of these special cases 
will immediately restore the problem of missing information. 
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If the matter on the brane has a perfect-fluid energy-momentum tensor, the local conservation equations 
([39]) reduce to 

p + <d(p + P ) = 0, 
D fl p+(p + p)A ll = 0, 

while the nonlocal conservation equations ([|2]) and (^) reduce to 



U 



ian L in »j in a , D^Pu^ + A V V. 



and 

(44) 
(45) 

(46) 
(47) 



Equation (47) shows that |^] if — and the brane energy-momentum tensor has perfect fluid form, then the 
density p must be homogeneous. The converse does not hold, i.e., homogeneous density does not in general imply 
vanishing S^. A simple example is the Friedmann case: Eq. (E7n is trivially satisfied, while Eq. (46) becomes 



U + AHU = . 



(48) 



This equation has the 'dark radiation' solution 



U = U [ — 

a 



(49) 



If B^y = 0, then the field equations on the brane form a closed system. Thus for perfect fluid branes with homogeneous 
density and £^ v = 0, the brane field equations form a consistent closed system. However, there is no guarantee that 
the resulting brane metric can be embedded in a regular bulk. 

It also follows as a corollary that inhomogeneous density requires nonzero S^. For example, stellar solutions on the 
brane necessarily have £^ v 7^ in the stellar interior if it is non-uniform. Perturbed Friedmann models on the brane 
also must have £ M „ 7^ 0. Thus a nonzero £^ v is inevitable in realistic astrophysical and cosmological models. 

For a perfect fluid at very high energies, i.e., p> A, and for which we can neglect U (e.g., in an inflating cosmology), 
Eqs. (§|) and (§4|) show that (16 1 



w 



2\tot 



2w + l, 



ct + w + 1 



(50) 
(51) 



where w — p/p and Cg = p/p. Thus at very high energies on the brane, the effective equation of state and sound 
speed are stiffened. This can have important consequences in the early universe and during gravitational collapse. 



For example, in a very high-energy radiation era, w 
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the effective cosmological equation of state is ultra-stiff: 



w tot rj |, i n late-stage gravitational collapse of pressureless matter, w = 0, the effective equation of state is stiff, 
w tot an( j £ ne e ff ec tive pressure is nonzero and dynamically important. 



VI. GRAVITATIONAL COLLAPSE TO A BLACK HOLE ON THE BRANE 



The dynamics of gravitational collapse on the brane is not yet properly understood, given the complications that 
are introduced by local high-energy effects and nonlocal effects. In general relativity, it is known that nonrotating 
collapse to a black hole leads uniquely to the Schwarzschild metric, with no 'hair', i.e., no trace of the dynamics of 
the collapse process. On the brane, this is no longer true. The Schwarzschild metric cannot describe the final state 
of collapse 0,^|[l)- 

since it does not incorporate the 5-dimensional behaviour of the gravitational potential in the 
strong-field regime (the metric is incompatible with massive Kaluza-Klein modes). The appropriate metric on the 
brane to describe the black hole final state is not known, but a non-perturbative exterior solution should have nonzero 
Eftv in order to be compatible with massive Kaluza-Klein modes in the strong-field regime. In the end-state of collapse, 
we expect a nonlocal field which goes to zero at large distances, recovering the Schwarzschild weak-field limit, but 
which grows at short range. Furthermore, £^ v may carry a Weyl 'fossil' record of the collapse process. 

The known black hole solutions on the brane include the Schwarzschild solution, for which = 0, but also other 
solutions with nonzero nonlocal effects. A vacuum on the brane satisfies the field equations 
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Rp, = -£ MV , #V = = £%, V^ = 0. (52) 

It follows that Einstein- Maxwell solutions in general relativity produce vacuum solutions on the brane, where the elec- 
tromagnetic field is replaced by a nonlocal Weyl energy-momentum tensor UM- Examples are the Reissner-Nordstrom- 



like solution discussed below, and a Vaidya-like solution |19[ . The equations ( |52| ) form a closed system on the brane 
in the stationary case, including the static spherical case, for which 

6 = = w„ = a„ v , U = 0=Q fl = V fl „. (53) 

The nonlocal conservation equations reduce to 

±D M W + \UA^ + D"7V + A v V liV = , (54) 
and the general solution of this and the remaining brane field equations, with metric of the form, 

ds 2 = -F{r)dt 2 + F~ 1 (r)dr 2 + r W , (55) 



is 



then 



F = .-OM i (56) 



\ M pJ r \M p 2 J r 
£ l iv = -\ JL- ) \ [u^u v - 2r^r v + h^} , (57) 

where is a unit radial vector. 

This solution has the form of the general relativity Reissner-Nordstrom solution, but there is no electric field on 
the brane. Instead, the nonlocal Coulomb effects imprinted by the bulk Weyl tensor have induced a 'tidal' charge 
parameter q. This parameter depends on the mass M on the brane, i.e., q = q(M), since this is the source of bulk 
Weyl field (leaving aside the more complicated case where there may be additional Weyl sources in the bulk). In order 
to preserve the spacelike nature of the singularity, we need q < 0. This is in accord with the intuitive idea that the 
tidal charge strengthens the gravitational field, since it arises from the source mass M on the brane. By contrast, in 
the Reissner-Nordstrom solution of general relativity, q = +Q 2 weakens the gravitational field. Negative tidal charge, 
q < 0, means that there is one horizon on the brane, outside the Schwarzschild horizon: 



M 
M 2 " 

±v± p 




M 4 

-i p 

1-9— - 



M 2 M2 



> r s . (58) 



The tidal-charge black hole metric does not satisfy the far-field r~ 3 correction to the gravitational poten- 
tial [p p^j20| , p^ | , as in Eq. (|^), and therefore cannot describe the end-state of collapse. However, Eq. (|6|) shows 
the correct 5-dimcnsional behaviour of the potential (oc r~ 2 ) at short distances, so that the tidal-charge metric should 
be a good approximation in the strong-field regime. This is the regime where the Schwarzschild metric on the brane 
fails, whereas the Schwarzschild metric is a good approximation at large distances. 

The bulk solution corresponding to the tidal-charge brane metric is not known, although numerical investigations 
have been performed pl|] . For the Schwarzschild special case q = 0, the brane is the surface z = 1 in AdSs, in 
conformal coordinates ]17|: 

d * 2 = {ix)it ds2+dz2 ^ (59) 

Here ds 2 is the 4-dimensional Schwarzschild metric, given by Eqs. (|5^) and (|56|) with q — 0. The singularity r — is 
a line along the z-axis, so that this bulk metric describes a 'black string'. The 5-dimensional horizon is the surface 
gtt = in the bulk, which is r = 2M /M 2 . This bulk horizon is a sphere of radius r = r s on each z =constant surface, 
so that it has a cylindrical shape in the z-direction. The bulk is singular at the AdSs horizon, and the black string 
horizon is unstable Jlj],^2| . 

If the physically realistic black hole bulk metric has the form 

ds 2 = -F{r, y)dt 2 + J(r, y)dr 2 + L{r, y)r 2 dQ 2 + dy 2 , (60) 

in Gaussian normal coordinates, then the bulk horizon is the surface F(r, y) — 0. Perturbative studies |]l8|,^2) and 
exact lower-dimensional solutions p3| have been used to probe the shape of the bulk horizon, but a non-perturbative 
5-D solution has not yet been found. 

Clearly the black hole solution, and the collapse process that leads to it, have a far richer structure in the brane- world 
than in general relativity. Another intriguing issue is how Hawking radiation is modified by bulk effects pj]. 
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VII. EXACT COSMOLOGICAL SOLUTIONS 



By Eqs. (f49|) and (pi]), the generalized Friedmann equation on a spatially homogeneous and isotropic brane is p5 26 



1 „2, 



re 2 A \ a 



(61) 



where if = 0, ±1. The nonlocal term that arises from bulk Coulomb effects (also called the dark radiation term) is 
strongly limited by nucleosynthesis [p5|j27|] : 



1 fU_ 

H 2 X \ p 



< 0.005. 



(62) 



cl 



This means that the nonlocal term is always sub-dominant during the radiation era, and rapidly becomes negligible 
after the radiation era. 

If we neglect the nonlocal term, and take A = = K , then the Friedmann equation during very high-energy inflation 
and early radiation-domination becomes 



k 2 n 2 
Tf 2 ~ P 



(63) 



The enhanced expansion rate, compared to general relativity (where H oc ^/p), introduces important changes to 
the dynamics of the early universe J^,^8|,^9| and leads to an increase in the amplitude of scalar || and tensor ]3C| l 
perturbations generated by inflation. If p 3> A at the start of the radiation era, then the solution of the Friedmann 
equation gives a oc i 1 / 4 . 

The generalized Raychaudhuri equation (|7(]) (with A = = U) reduces to 



H + H 



6 ^ 



p + 3 P +(2p + 3 P )^ 



and shows that the condition for inflation on the brane (a > 0, i.e. H + H 2 > 0) is || 

1 f2p + X" 



w < 



p + X 



(64) 



(65) 



As p/X — > 0, we recover the general relativity result, w < —h. But in the very high-energy regime, we find w < 



The dynamics of homogeneous but anisotropic Bianchi branes has also been investigated |14j,|29j]. In particular, 
high-energy bulk effects can strongly alter the dynamics near the singularity, because matter can dominate over shear 
anisotropy, opposite to the case of general relativity fl4j| . 

The bulk metric for a flat Friedmann brane may be given explicitly pq| : 



dS 2 = -N(t, yfdt 2 + A(t, yfdx 2 + dy 2 , 
where A(t, 0) = a{t) and t is proper time on the brane, so that N(t, 0) = 1. The metric functions are 



N 



A 
a 



A 2 = —^ a 2 - 
2A V X) 

2p 



m n at 



.I2 + - + -U 



X 2 a 2 



cosh(2^y) 



(66) 



(67) 



(68) 



with p, = ^^/Q/X. This bulk metric is Schwarzschild-AdSs spacetime [[L3| (in Gaussian normal coordinates), with 
the mass parameter of the black hole in the bulk proportional to IA . 
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VIII. PROPAGATION AND CONSTRAINT EQUATIONS 



Section V gave propagation equations for the local and nonlocal energy density U and flux Q M . The remaining 
covariant equations on the brane are the propagation and constraint equations for the kinematic quantities O, A lll 
Cjiii/j an d for the nonlocal gravitational field on the brane. The kinematic quantities govern the relative motion 
of neighbouring fundamental world-lines. The nonlocal gravitational field on the brane is given by the brane Weyl 
tensor C^ va p. This splits into the gravito-electric and gravito-magnetic fields on the brane: 

Efiv — C t _ iav pu a u 13 — , H^y = \e ^apC^ v^u 1 = , (69) 

where EL„ must not be confused with The Ricci identity for and the Bianchi identities V^C^q,^ — 

V[ M (— R u \ a + \Rg v ]ci) produce the fundamental evolution and constraint equations governing the above covariant 
quantities [ gj] . The field equations are incorporated via the algebraic replacement of the Ricci tensor i? M „ by the 
effective total energy- momentum tensor, according to Eq. (|3l|). The brane equations are derived directly from the 
standard general relativity versions [p2"| by simply replacing the energy- momentum tensor terms p, ... by p tot , .... For 
a perfect fluid or minimally-coupled scalar field, the general equations Q reduce to the following. 

Generalized Raychaudhuri equation (expansion propagation): 

6 + |e 2 + a^a^ - 2lu^ - D^A^ + + \n\p + 3p) - A 

= -^(2p + 3p)^---^U. (70) 



Vorticity propagation: 



j {fi) + |6w M + fcurl^ - a^Lu" = . (71) 



Shear propagation: 
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(itv) + + E^ u - ~D^A V ) + <7 a ^(T v ) a + uj^Wv) - A^A V ) = -^-V^v . (72) 



Gravito-electric propagation (Maxwell- Weyl E-dot equation): 

E(n V ) + QEfj, u - curXH^ + \n 2 {p + p)a liV 
- 2A a £ aj} ^H v ^ - 3<J a ^E v ) a + uj a e al3 ^E I/ ^ 

A K A L 

+ &A (li Q v) + 3a a {fl T„ )a + 3uj a s a ^V^] . (73) 
Gravito-magnetic propagation (Maxwell- Weyl H-dot equation): 

H (llv) + QH„ V + cuxlE^ - 3a a{ „H u) a + ^e^H^ + 2A a e a ^E u / 
= -|- [curlTV - 3lu^Q„ } + ^^e^pQ 13 ] . (74) 



K 2 A 



Vorticity constraint: 



Shear constraint: 



- A^uj^ = . (75) 



D"v ~ cur1 ^ - ! d m© + 2[w, A}„ = — l~Q„ . (76) 

Gravito-magnetic constraint: 

curler^ + D^u v ) - + 2A( fl u> u ) = . (77) 
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Gravito-electric divergence (Maxwell- Weyl div-E equation) 



TTEp, - i« 2 D M(0 - [a, H]„ + 3H^u v 

K 2 P 1 



D mP + — (2D M W - 2QQ, t - 3D"?V + Zct^Q" - 9[w, Q] M ), (78) 



where the covariant tensor commutator is [W, Z]^ — e lia pW a . y Z^' y . 
Gravito-magnetic divergence (Maxwell- Weyl div-H equation): 

D V H„ V ~k 2 { P + P )uj^ + [a, E]„ - ZE^uj v 



= ^(p + p)^ + (8Ww„ - 3curl Q M - 3[a,^ t - 3^0;") . (79) 



Gauss-Codazzi equations on the brane (oj„ = 0) JT3]: 



#<U) + °"<M"> + 0O > ~ D <M j4 f> ~ ^(m^> = "TT^V > ( 80 ) 
+ |9 2 - <7 M ^" - 2k 2 p - 2A = ^ + (81) 

J A K Z X 



where Rh v is the Ricci tensor for 3-surfaces orthogonal to u M on the brane and i? x = WR 



The standard 4-dimcnsional general relativity results are regained when A -1 — * 0, which sets all right hand sides 
to zero in Eqs. (|70|)— (|8~l]) . Together with Eqs. (flij)-(|l7|), these equations govern the dynamics of the matter and 
gravitational fields on the brane, incorporating both the local (quadratic energy-momentum) and nonlocal (projected 
5-D Weyl) effects from the bulk. Local terms are proportional to p/X, and are significant only at high energies. 
Nonlocal terms contain U, Q M and V^v, with the latter two quantities introducing imperfect fluid effects, even though 
the matter has perfect fluid form. 

Bulk effects give rise to important new driving and source terms in the propagation and constraint equations. The 
vorticity propagation and constraint, and the gravito-magnetic constraint have no direct bulk effects, but all other 
equations do. Local and nonlocal energy density are driving terms in the expansion propagation. The spatial gradients 
of local and nonlocal energy density provide sources for the gravito-electric field. The nonlocal anisotropic stress is a 
driving term in the propagation of shear and the gravito-electric/ -magnetic fields, and the nonlocal energy flux is a 
source for shear and the gravito-magnetic field. The Maxwell- Weyl equations show in detail the contribution to the 
nonlocal gravito-electromagnetic field on the brane, i.e., (E^, H^), from the nonlocal 5-dimensional Weyl field in 
the bulk. 



The system of propagation and constraint equations, i.e. Eqs. (44)-(|47j) and (|70[)-(|8l|), is exact and nonlinear, 
applicable to both cosmological and astrophysical modelling, including strong-gravity effects. In the next section we 
will linearize the system in order to study cosmological perturbations on the brane. A different linearization scheme 
could be developed for studying compact objects. 

In general the system of equations is not closed: there is no evolution equation for the nonlocal anisotropic stress 
V^v As noted above, closure on the brane can arise in special cases, when the nonlocal conservation equations 
and (filil ) are sufficient to determine E^ v . These special cases include: 

Vfn> = 0, as on a Friedmann brane; 

&iiv = and D^p = 0, with perfect fluid matter, as in special cases of Friedmann and Bianchi branes; 
stationary brane, with £ M „ = 0, so that Eqs. ( [42| ) and ( f43| ) reduce to constraint equations; 
branes with isotropic 3- Ricci curvature, ui^ = = Rf^ v ) ■ 

The last case follows from Eq. d8C|), which becomes an equation defining V^v in terms of quantities already de- 
termined via other equations. Special cases are Friedmann and Bianchi I branes. In general, for anisotropic 3-Ricci 
curvature, Eq. j8Q ) is an equation determining Rf^ in terms of other quantities. 

Although the special cases above can give consistent closure on the brane, there is no guarantee that the brane is 
embeddable in a regular bulk. This is the case for a Friedmann brane, whose symmetries imply (together with Zi 
matching) that the bulk is Schwarzschild-AdSs ]l3[| . A Schwarzschild brane can be embedded in a 'black string' bulk 
metric, but this has singularities |l7|]2i||. 
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IX. COSMOLOGICAL PERTURBATIONS 



The dynamics of the background homogeneous brane provides important constraints on brane-world cosmologies. 
Roughly speaking, these constraints amount to the statement that the brane-world should reproduce the standard 
Friedmann dynamics from nucleosynthesis onwards. However, much stricter constraints are implied by the growing 
body of data on the cosmic microwave background (CMB) radiation and large-scale structure (LSS), since this data 
indirectly probes earlier times. In order to confront brane-world models with the data, we need to study perturbations 
on the brane. 

As is apparent from the exact nonlinear equations in the previous section, cosmological perturbations on the brane 
are qualitatively more complicated than in general relativity. Not only is the background more complicated, but 
effects from the bulk that are imprinted on the brane include degrees of freedom that cannot be predicted from 
data available to brane observers. A complete understanding of brane perturbations therefore necessarily involves 
the analysis of bulk perturbations - for which the mode equations are partial differential equations, with nontrivial 
boundary conditions. 

A covariant brane-perturbation theory has been developed Q| and applied to large-scale density perturbations [fl6| . 
Metric-based general formalisms for bulk perturbations have been developed |33| , and large-scale perturbations gener- 
ated from quantum fluctuations during de Sitter inflation on the brane have also bee 
Large-scale scalar perturbations and their impact on the CMB have been analysed 

The results so far may be summarized as follows: 



i partly computed |) 34 3£ 
ffi\. 



Tensor perturbations: Bulk effects produce a massless mode during inflation and a continuum of massive Kaluza-Klein 
modes HQ, with m > \H . The massive modes stay in the vacuum state, and on large scales there is a constant 
mode with enhanced amplitude (compared to general relativity) J3(| . We expect no qualitative change on large scales 
to 4-dimensional general relativity tensor modes in the CMB and LSS, but there could be a significant change on 
small scales due to the massive modes. 

Vector perturbations: Bulk effects can support vector perturbations, even without matter vorticity P,p6|, but these 
modes, which do not arise in general relativity, are massive (there is no normalizable massless modejTand stay in 
the vacuum state during inflation on the brane [ p6| . The momentum can be determined on large scales without 
solving the bulk perturbations, but the vector Sachs-Wolfe effect cannot be found on-brane p6[ , because of the non- 
local anisotropic stress, which is undetermined on the brane. There are possibly qualitative changes to 4-dimensional 
general relativity vector modes in the CMB and LSS on large scales. 

Scalar perturbations: Bulk effects introduce a non-adiabatic mode on large scales P,|l6| , p7[ . Density perturbations 
on large scales can be solved on-brane, without solving for the bulk perturbations llqH but the Sachs- Wolfe effect 
cannot be found on-brane [ p7[ , because of the nonlocal anisotropic stress, which is undetermined on the brane. There 
are possibly qualitative changes to 4-dimensional general relativity scalar modes in the CMB and LSS, even on large 
scales, and probably significant changes on small scales. 

The quantitative results up to now are confined to large scales. Further progress requires the integration of the 
bulk mode equations. In the simplest case, for tensor perturbations on a flat Friedmann brane, the mode equation 
is @ 

where A and N are given by Eqs. (^) and (|68|). The perturbed metric, in Gaussian normal coordinates, is 

ds 2 = -N(t, yfdt 2 + A(t, y) 2 [<S y - + U tj {t, x, y)] dx i dx J + dy 2 , (83) 
where Hi* = = d^Hij and 

(t, x, y) = H(t, y) exp(i£; • x) e tj , (84) 
with Cij a polarization tensor. At the brane, the boundary condition is 



m 

dy 



(85) 

v=o 
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where Tr T (t) exp(ifc • x) eij is the anisotropic stress of matter on the brane. 

In order to illustrate some of the features of perturbations on the brane, we can use the brane-based covariant 
approach Q]. The exact nonlinear equations in previous sections can be linearized as follows. The limiting case of the 
background Fricdmann brane is characterized by the vanishing of all inhomogeneous and anisotropic quantities: 

= V^ = W^=0, (86) 

where / = and — A^^u)^, Q^, and W^u — cr^, E^, H^jV^. These quantities are then first-order of 

smallness in the linearization scheme, and since they vanish in the background, they are gauge-invariant B^j . The 
linearized conservation equations (assuming adiabatic matter perturbations) are 

p + e(p + P ) = 0, (87) 

^+(^^, = 0, (88) 

U + §6W + D' 1 Q / , = ) (89) 

Qfi + 4H Q M + ±D M W + \UA^ + YfV^v = -\n\p + p)D M p . (90) 

Linearization of the propagation and constraint equations leads to: 



6 

A k~\~ 

n + 2Hu M + icurl A„ = , (92) 



e + ±e 2 -jy^ + ± K 2 (p + 3 P ) - a = -i K 2 (2 P + 3p)^ - — u, (9i) 



+ 2Ha^ + E^ - B (fi A u) = -^TV , (93) 



ft 2 A 

E^ v + iHE^ v - curl-ff^ + \n 2 {p + p)<7^ v = -\n 2 {p + p)^<V 



' AU(T^ + 2P( Jlu) +&P ia , + 2D {ll Q v) , (94) 



K 2 X 

+ 3HH^ + curlE^ = -4r curl7 V . ( 95 ) 



D"^ = 


0, 




DV^ - 


- curlw M - fD^e = - 


— Q 

K 2 A 


curlcr Ml , 


+ D^w^) — if M „ = , 






1 2 n k2 Pt, 


1 

f ^A 


D"^ - 


- K 2 (p+p)uJ^ = K 2 (p 





(96) 
(97) 



[2D^W - 20e M - SD"^] , (99) 

^ + 4t [8Ww m - 3curl Qfj] . (100) 

Equations (^), (89) and (^lj) do not provide gauge-invariant equations for perturbed quantities, but their spatial 
gradients do. 

A covariant, gauge- invariant (and purely local) splitting into scalar, vector and tensor modes is given by (com- 
pare@) 

V» = D M V + % , (101) 
WV - D^Dy) W + B {li W v) + , (102) 

where an overbar denotes a transverse (divergence- free) quantity (note that W^ v is already tracefree). Purely scalar 
modes are characterized by 

% = % = %u = , (103) 

and standard identities [Q, the vorticity constraint equation ( p6| ) and the gravito-magnetic constraint equation ( pSj ) 
then show that 

curl Vp = = curl W„ v , D V W^ = §D 2 (D^W) , u„ = = . (104) 
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Vector modes obey 

V„ = %, = T> (fl W v) , curlD M / = -2/u; M , 

and it then follows that 

V V W^ = ±D 2 W> M , curlW^ - ^ {IX W V) . 
Tensor modes are covariantly characterized by 

D M / = = V M) W IIV = W IU> . 



(105) 
(106) 
(107) 



A. Density perturbations on the brane 

We define the density and expansion scalars (as in general relativity J37J) 

A = —D 2 p, Z = a 2 V 2 Q, 
P 

and scalars describing inhomogeneity in the nonlocal quantities: 



(108) 



U=— D 2 U, Q = -D 2 Q, P=-D 2 V. 



(109) 



Then we take the comoving spatial Laplacian of Eqs. (j87|), ( p9[) and (|9l|), using standard identities H). This leads to 
a system of 4 evolution equations for A, Z, U, and Q. In general the system is under-determined, since there is no 
evolution equation for P. However, P only arises via its Laplacian, and on large scales, the system closes and brane 
observers can predict density perturbations from initial conditions intrinsic to the brane. The system on large scales 
is 



A = 3wHA - (l + w)Z 
Z = -2HZ 



6 P tt 12 



. p Ylc 2 s U 
1 + 4 + 3tu)£ - s -—r- 



A. 



U = (3w-l)HU- 



1 + w 



(110) 
(111) 

(112) 

(113) 

In standard general relativity, only the first two equations apply, with A" 1 set to zero in Eq. (Ill), so that we can 
decouple the density perturbations via a second-order equation for A, whose independent solutions are adiabatic 
growing and decaying modes. Bulk effects introduce a new non-adiabatic mode: there are 3 coupled equations in A, Z 
and U, plus a decoupled equation for Q, which is determined once the other 3 quantities are solved for. Some solutions 
are found in |l6| ] for U = in the background, and they show how the fluctuations U in U introduce effective entropy 
perturbations, and how A evolves differently from the standard general relativity case when p^> \ (see also p7j]). 



= (l-3w)HQ-^-U+^- 
3a ba 



S>C S U 4 2/1, \ 

— r K p (1 + W) 



A. 



B. Vector perturbations on the brane 



The linearized vorticity propagation equation (92) does not carry any bulk effects, and vorticity decays with 
expansion as in standard general relativity, reflec ting the fact that angular momentum conservation holds on the 
brane. The gravito-magnetic divergence equation (|100|), becomes 



D 2 i^ = 2k 2 ( P + p) 



1 



K 2 A 



[8Ww M -3curlQ M ] , 



(114) 



on using Eq. ( JlOq ). This shows that the local and nonlocal energy density and the nonlocal energy flux provide 
additional sources for the brane gravito-magnetic field. In standard general relativity, it is necessary to increase the 
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angular momentum density (p + p)^ in order to increase the gravito-magnetic field, but bulk effects allow additional 
sources. In particular, unlike in general relativity, it is possible to source vector perturbations even when the vorticity 
vanishes, since curl may be nonzero. 

W e can find a closed system of equations on the brane for uj^ and curlQ M , on large scales. Using Eqs. ([9Cj), 
f05|) and fll06|), we find 



K 



l-K) Ha^=0, 
(1 - Zw)HP M = 



U 



(115) 
(116) 



where 



auj b 



0M 



a i n 
-curl Q, 

P 



(117) 



are dimensionless vector perturbations. Thus we can solve for the matter and nonlocal vorticity on large scales, 
without solving the bulk perturbation equations, similar to the case of scalar perturbations. The matter vorticity 
sources nonlocal vorticity, but during slow-roll inflation this only leads to a highly damped mode oc a~ 8 . In very 
high-energy radiation domination however, there is a weakly growing mode: 



Bb^ In I — 



(118) 



where B 



C. Tensor perturbations on the brane 



For the transverse traceless modes on the brane, the scalar equations reduce to background equations, the vector 
equations drop out and we can derive a covariant wave equation for the tensor shear from the remaining equations j9^ ) 
© and (|§): 
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2A + W(p-Zp) 



2HV 



[US 



(119) 



In standard general relativity, the right hand side falls away, and once a^ v is determined, one can determine 
from Eq. ( p3| ) and H^ v from Eq. (p8| ) (the shear is a gravito-electromagnetic potential) . Nonlocal bulk effects provide 
driving terms that are like anisotropic stress terms in general relativity p9[ . In the latter case however, the evolution 
of anisotropic stress is determined by the Boltzmann equation or other intrinsic physics. The nonlocal anisotropic 
stress Vfn, from the bulk Weyl tensor is not determined by brane equations, and can in principle introduce significant 
changes to the shear via Eq. (119). 



X. CONCLUSION 



By adopting a covariant approach based on physical and geometrical quantities that are in principle measurable by 
brane-observers, we have given a classical analysis of intrinsic cosmological and astrophysical dynamics in generalized 
Randall-Sundrum-type brane-worlds. We have emphasised the role of the Weyl tensors in the bulk and on the 
brane, and we have carefully delineated what can and can not be predicted by brane observers without additional 
information from the unobservable bulk. Intriguing new effects are introduced by the extra dimension, and there are 
exciting problems to be solved, especially on branc-world black holes and gravitational collapse, and on cosmological 
perturbations and the CMB and LSS. 
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